Constraining Isocurvature Initial Conditions with WMAP 3-year data 



o 
o 

(N 



Rachel Bean', Joanna Dunkley^''', Elena Pierpaoli ^ 
" Dept. of Astronomy, Cornell University, Ithaca, NY, USA. 
t Dept. of Physics, Princeton University, Princeton, NJ, USA. 
^ Dept. of Astrophysical Sciences, Princeton University, Princeton, NJ, USA. 
* Physics and Astronomy Department, University of Southern California, Los Angeles, CA 90089-04-84, 

We present constraints on the presence of isocurvature modes from the temperature and polariza- 
tion CMB spectrum data from the WMAP satellite alone, and in combination with other datasets 
including SDSS galaxy survey and SNLS supernovae. We find that the inclusion of polarization data 
allows the WMAP data alone, as well as in combination with complementary observations, to place 
improved limits on the contribution of CDM and neutrino density isocurvature components individ- 
ually. With general correlations, the upper limits on these sub-dominant isocurvature components 
are reduced to ~60% of the first year WMAP results, with specific limits depending on the type 
of fluctuations. If multiple isocurvature components are allowed, however, we find that the data 
still allow a majority of the initial power to come from isocurvature modes. As well as providing 
general constraints we also consider their interpretation in light of specific theoretical models like 
the curvaton and double inflation. 
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I. INTRODUCTION 

Impressive recent developments in measurements of 
the Cosmic Microwave Background (CMB) temperature 
and polarization anisotropies P, 0, |3|, large scale 
structure (LSS) [1, H, 0, H, supernovae [13, IHI and 
the Lyman-Qf forest have enabled stringent test- 
ing of the cosmological model including the composition 
of the universe and the initial conditions that seeded 
inhomogeneities. Although the simplest initial condi- 
tions, arising from single field inflation, predict scale- 
invariant, adiabatic inhomogeneities, this is far from the 
only possibility. Isocurvature modes are predicted by 
a wide r ang e of scenarios, for example multi-field infla- 
tion [H, [11 [13, topological defects [ll[ll, and 
through the decay of particl es p rior to nucleosynthesis 
such as a scalar curvaton [2l|, [H, [H, [H, [2^, [2^ , ax- 
ions [13] or the Affleck-Dine model of baryogenesis [l3] ■ 
The degree of correlation with adiabatic perturbations 
can vary across the full spectrum, from completely (anti) 
correlated in the curvaton scenario, intermediate correla- 
tion in some multi-field inflation models, to uncorrelated 
modes in cosmic string scenarios. Most proposed scenar- 
ios generate solely baryon or Cold Dark Matter (CDM) 
isocurvature modes [28| ; mechanisms generating neu- 
trino density and velocity isocurvature modes are also 
possible [13, [lOl , though the latter are more difficult to 
motivate. 

Although the simplest adiabatic scenario is in com- 
plete agreement with the current data the question 
still arises, how large a contribution could isocurvature 
modes make and how sensitive is the cosmological param- 
eter estimation to their inclusion? In this paper we con- 
front these questions, specifically focusing on constraints 
from the WMAP CMB temperature and polarization 
power spectra [H, [1 i, [1, SDSS galaxy matterp ower 
spectrum [^ and the SNLS supernovae survey Ex- 
periments in the past ten years have ruled out models 
with purely isocurvature perturbations [3l|, [H, [H, [s^] , 



but have not excluded those with admixtures of adiabatic 
and isocurvature modes. Constraints placed on models 
with a single isocurvature mode, in light of the first year 
WMAP results, indicate that it must be sub-dominant 
[H, HE [13, H [H, HE [U, m, but models with multi- 
ple modes may have larger non-adiabatic contributions 
pol [43l . . Here we analyze the constraints in light of 
the 3-year WMAP data release of temperature and polar- 
ization data for which only an analysis for perfectly cor- 
related CDM or baryon isocurvature perturbations has 
been conducted so far (l^ . lisj . 

In section |ll] we outline the approach and parameter- 
izations used in the analysis. In sections IIIII and IIVI we 
respectively establish the constraints on scenarios with 
purely correlated and generally correlated single isocur- 
vature modes. In doing so we refiect on what the limits 
on isocurvature contributions imply for some key par- 
ticle based theories which depart from purely adiabatic 
perturbations. In section |V] we expand our analysis to 
allow multiple isocurvature modes with general correla- 
tions, in which destructive cancellations can allow a large 
fractional isocurvature contribution. We conclude with 
a summary of our findings and implications for future 
cosmological observations in section IVIl 



II. APPROACH 

Perturbations to the metric may give rise to both cur- 
vature perturbations on comoving hypersurfaces, as well 
as entropy perturbations where the space-time curva- 
ture vanishes at early times. The former are termed 
adiabatic perturbations and may be quantified by the 
curvature perturbation, TZ. The latter are isocurvature 
modes, quantified by the entropy perturbation Sx ~ 
Spx/ {px +Px) — Sp-y/ {p^ +P-y) in the case of density per- 
turbations, Sp, between photons and a fiuid x, which may 
be CDM or baryons. There are two further isocurvature 
modes where the sum of the neutrino and photon den- 
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sities, or momentum densities, are initially unperturbed, 
whose initial conditions are given in [sol . [48| . 

To first order, initial conditions 5pc = —Spt, 6p-y = 
6pi, ~ allow a time independent solution to the per- 
turbation equations for the pressureless matter compo- 
nents. Baryon and CDM isocurvature initial conditions, 
up to a factor of flc/^b, are essentially observationally 
indistinguishable and therefore we only consider CDM 
isocurvature scenarios of the two here. 

The perturbations may be characterized using the co- 
variance matrix A, where 

\,ik)sHk~k') = (^Ay (;t,(k)x;(k')> (1) 

as shown in [soj . The subscript i may take the values 
{A,C,N,V}, labehng the adiabatic (AD), CDM isocur- 
vature (CI), neutrino isocurvature density (NID), and 
neutrino isocurvature velocity (NIV) components respec- 
tively, and the random variable Xi corresponds to the 
amplitude of the ith mode. 

We assume that the elements of the matrix A can be 
parameterized as 

Ay(fc)=Aj f^j , n,j^-{n,+nj) (2) 

in terms of a set of amplitude parameters Aij , power-law 
spectral indices rii, and the pivot value ko — 0.05 /Mpc. 

The angular power spectrum of a given correlation ij 
is given by 

Cr= ^Al{k)ei{k)&,{k), (3) 

where Oi{k) is the photon transfer function for initial 
condition i. An admixture of the adiabatic mode with 
a single isocurvature mode can be expressed in terms of 
the pure adiabatic, isocurvature and wholly correlated 
spectra, such that 

C, - AAACf^ + AiiC',' + 2AAiCf', (4) 

where / labels C, iV or V^. This is commonly parameter- 
ized using one of the following: 

Ce - Cf^ + B^C^' + 2B cos 9Cf^, (5) 
= (1 - a)C/^ + aCP + 2(3^/a{l - a)C/^ (6) 

to within an overall normalization factor, with the former 
used by 'si', 3?, "36^ and the latter by [33, S . Using 
equation ([5]) , the amplitude and correlation phase of the 
isocurvature contribution are given by a and /? respec- 
tively. These in turn can be related to the overall ratio of 
isocurvature to adiabatic component given by B = S/TZ, 
with a = B'^ /{l + B"^), and general correlation /3 — cosd. 
We will use a and P as parameters in this analysis, for 
models with a single isocurvature mode, sampling with 
f3 directly rather than 2/3-^ a(l — a). With these defini- 
tions, a positive correlation between the adiabatic and 



isocurvature perturbations will produce a wholly corre- 
lated CMB spectrum with a negative amplitude at large 
scales, in agreement with e.g. [sl, [SI, [H, [1^ and the 
CAMB package. The primordial adiabatic perturbation 
may be defined such that these correlated spectra have 
the opposite sign, as used in [13, |4l|, [H, . 

The above parameterizations do not naturally extend 
to the addition of more than one isocurvature mode, for 
which a method is given in [i^ 0| . Here 

N 
ij = l 

for N perturbation modes, where ^fj^i zf^ — 1, and 
any matrix z with a negative eigenvalue is assigned a zero 
prior probability. Cg are normalized to have equal CMB 
power in each mode i, such that quan tify the physi- 
cally observable power in each mode^ [55|. By definition 
the coordinates lie on the surface of a A'^^-dimensional 
unit sphere, which can be sampled uniformly with N'^ — 1 
amplitude parameters using the volume preserving map- 
ping shown in [40| . To quantify the isocurvature contri- 
bution in the case of multiple modes, a measure is given 
by Tiso = Ziso/{ziso + zaa), where Zi^o = ^/i^^^^^- 
This is equivalent to the fiso parameter in l40|. l43l . |44| . 
but different to the fiso used in e.g. [H, l36l. l42l| . which 
corresponds to B in equation ([5]). 

It is important to note that constraints on a can de- 
pend strongly on the pivot scale, whereas the Zij pa- 
rameters are independent of the pivot. As such, the Zij 
parameterization is a useful direct measure of isocurva- 
ture, whereas a can be misleading when the isocurvature 
and adiabatic spectral tilts are allowed to vary indepen- 
dently, as will be discussed in section IIVI For ease of 
comparison with previous analyses, however, for single 
components we present constraints for {a, /?} (pivoted at 
k = 0.05/Mpc) and give equivalent constraints on {%}, 
since the relation between the single mode parameteri- 
zation {a,P} and the {zij} is not trivial. For multiple 
modes we solely present results using the {zij} parame- 
terization. 

We parameterize our cosmological model in terms 
of a ACDM scenario using the following parameters: 
flbh^, ftch^, T, J7a, bsDSSi and an overall scalar amplitude 
parameter, limiting our search to flat models with scalar 
fluctuations and assume 3.04 massless neutrinos species 
with zero chemical potential. We do not investigate here 
broader parameter spaces including tensor modes, run- 
ning in the scalar spectral index, spatial curvature or 
evolving dark energy. Throughout this paper we use a 
single parameter for both the isocurvature and adia- 
batic spectral index. However, we do relax this constraint 
in section HVl where we discuss the implications of a sin- 
gle CDM isocurvature mode generally correlated with the 
adiabatic one. 

We find constraints using the 3-year WMAP CMB 
data alone P, [^ for CDM and neutrino density isocur- 
vature scenarios with a variety of degrees of correla- 
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tion, and constraints on one, two and three isocurva- 
ture components more generally in combination with 
data from the SDSS galaxy survey , including a Gaus- 
sian prior on the SDSS bias measurement 1.03 ± 0.15 fi\ 
and non-linear corrections [46] , and supernova data from 
SNLS We include BBN estimates of the baryon 

to photon ratio, conservatively encompassing measure- 
ments of both Deuterium (ryio = 6.4 ±0.7) and Helium-3 
(?7io = 6.0 ± 1.7) by imposing a Gaussian prior of 
rife/i^ = 0.022 ± 0.006. Small-scale polarization data do 
not yet noticeably tighten the constraints so we do not 
include them in the analysis. As we discuss in the analy- 
sis and conclusion, however, future small scale polariza- 
tion data could well be an important test of isocurvature 
scenarios. We generate CMB and matter power spectra 
using the CAMB package [i^, which is consistent with 
our perturbation definitions. The likelihood surfaces are 
explored using Markov Chain Monte Carlo methods, ap- 
plying the spectral convergence test described in [4§| , and 
the Gelman and Rubin convergence test [s^] ■ We also use 
a downhill simplex method to find the best-fit likelihood, 
starting from the maximum likelihood point sampled by 
the chains, since the likelihood peak is only sparsely sam- 
pled by MCMC in high dimensional spaces. 



III. UNCORRELATED AND PERFECTLY 
CORRELATED ISOCURVATURE: SINGLE 
MODES 

In a variety of theoretical scenarios the isocurvature 
and adiabatic fluctuations in matter can arise out of 
a single mechanism and subsequently have well-defined 
degrees of correlation. We first consider constraints 
on purely correlated (/? = 1), uncorrelated (/? = 0) 
and anti-correlated (/? = —1) CDM and neutrino den- 
sity isocurvature fluctuations from WMAP data alone 
(TT-I-TE+EE), and in combination with other datasets. 
Figure [1] shows that for both CDM and neutrino den- 
sity isocurvature, purely correlated scenarios are the 
most tightly constrained by the data. For CDM modes 
with WMAP+SDSS-l-SNLS-hBBN (and WMAP only) 
we find a < 0.009(< 0.01) at 95% confidence limit 
(c.l.) for purely correlated and < 0.009(< 0.045) for 
anti-correlated modes (consistent with [12l |. who find 
a < 0.005 including Ly-a data). The corresponding lim- 
its for neutrino density modes are a < 0.017(< 0.025) 
and < 0.026(< 0.083). The wholly uncorrelated isocurva- 
ture modes are allowed to contribute a much more signif- 
icant fraction of the overall power, with a < 0.11(< 0.26) 
and a < 0.21 (< 0.55) (95% c.l.) for CDM and neutrino 
density isocurvature modes. We consider the statistical 
support for presence of isocurvature flucutations by us- 
ing the bestfit likelihood , £, calculated by comparing 
the theoretical spectrum predicted by the cosmological 
scenario, {C*'*}, to the observed temperature and polar- 
ization maps or spectra, x, in light of the statistical and 
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FIG. 1: 1-dimensional likelihood distributions for CDM (top) 
and neutrino (bottom) density isocurvature models with per- 
fectly correlated j3 = 1, anti-correlated /3 = — 1 and uncor- 
related /3 = isocurvature perturbations using WMAP only 
data (black full) and combined with SDSS, SNLS and BBN 
data (red dashed). 

systematic errors encoded in the covariance matrix C, 

/:(x|c-)= n ^^Pi^=^. (8) 

J 7 t vdet 6 

datasets 

The 3-year WMAP data is entirely consistent with no 
isocurvature contribution being required, having no im- 
provement over the bestfit likelihood for the adiabatic 
scenario — 21n£ = 11252 ( arising from the joint pixel/ 
spectrum based likelihood approach outlined in [l|, 01 
with total = E^C^^x = 3279 for 3244 degrees of 
freedom and E Indet C = 7973 ). 

For WMAP data alone significant cosmological pa- 
rameter degeneracies exist between a and ilc/i^ and ns, 
arising because the principal effects of the isocurvature 
modes are modifications to the large scale temperature 
fluctuations. We find no significant degeneracy between 
the isocurvature fraction and the optical depth to reion- 
ization. The impact of these degeneracies are most signif- 
icant for the uncorrelated and anti-correlated CDM and 
neutrino density modes while they have only a nominal 
effect on the correlated mode constraints. For the uncor- 
related and anti-correlated modes Slc/i^ is decreased by 
roughly 1 cr, and Us is increased by 2 cr from the adiabatic 
value. The addition of SDSS, SNla and BBN datasets 
tighten the constraints by truncating these degeneracies 
and bringing the values back towards the fiducial values. 
The correlated CDM and nuetrino modes are little im- 
proved by the inclusion of SDSS -f SNLS data because 
in these scenarios there is no significant degeneracy be- 
tween a and and ^Ich^, the two parameters that are 
significantly better measured by the inclusion of the com- 
plementary datasets to the CMB. 

These constraints have implications for the curvaton 
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scenario |24l |. which includes accelerated expansion by 
inflation but allows for primordial perturbations to be 
generated by the decay of a distinct scalar field, the cur- 
vaton. While no unique prescription for the generation 
of fluctuations in the curvaton scenario exists, there are 
a range of scenarios where the curvaton gives rise to the 
cold dark matter isocurvature perturbations, which in 
general predict Uadi = n-iso- The curvaton scenario does 
not provide a unique prescription for the generation of 
fluctuations, however in its simplest form it predicts the 
existence of cold dark matter isocurvature perturbations 
with Uadi = n.iso- This is because the curvature and en- 
tropy perturbations are related to the gauge invariant 
Bardeen variable f 



7^ = 



(9) 
(10) 



The kind and amount of isocurvature depends on when 
the curvaton field decays and CDM is created [5l|. Sce- 
narios in which CDM is generated prior to curvaton decay 
have S,cDM = and the entropy and adiabatic fluctua- 
tions perfectly anti-correlated {(3 = 1) yielding a = 0.9 
which remains ruled out at high significance. If the 
CDM is to be generated by the curvaton decay then 
[3 — ^1 and the amount of isocurvature refiects the ratio 
of the curvature fiuctuation after decay to before it, r = 
(,bef ore /Rafter which using the suddcn decay approxima- 
tion r « Pcurvaton/ Plot, r = [l + \^ [U . Our 

analysis therefore sets limits on the curvaton decay, with 
0.97 < r < 1 (95% c.l.) from WMAP+SDSS-hSNLS, 
comparable to the ones obtained by Beltran et. al. jl^l 
(r > 0.98 95% c.l.), who also included Lyman a con- 
straints and used a slightly different parametrization. For 
neutrino isocurvature modes generated by density pertur- 
bations we find r > 0.94 (95% c.l.). The most practical 
mechanism, however, for generating neutrino isocurva- 
ture perturbations is through a perturbation in the lepton 
number f23|,and subsequent non-zero chemical potential, 
not analyzed here. 



IV. GENERALLY CORRELATED 
ISOCURVATURE: SINGLE MODES 

We next consider constraints on generally correlated 
isocurvature fluctuations from combined WMAP, SDSS, 
SNla and BBN data, including the CDM density, neu- 
trino density and neutrino velocity modes individually. 

Constraints on the two-dimensional isocurvature am- 
plitude and correlation spaces are shown in Fig. [5] and 
summarized in Table HI For CDM isocurvature we find 
a < 0.15 at the 95% c.l. and no overall improvement 
in the goodness of fit — 21n£ = 11383. Neutrino den- 
sity models have a < 0.18, and neutrino velocity models 
a < 0.26. The CDM mode prefers a small positive cor- 
relation with the adiabatic mode. 
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TABLE L 95% upper (or lower) limits, or means and 68% 
confidence levels, for mode contributions for models with gen- 
erally correlated isocurvature. Scenarios using WMAP, SDSS 
and SNLS data in which the adiabatic and isocurvature spec- 
tral indexes are both fixed to be identical and where they are 
allowed to differ are shown. The best-fit likelihood and the 
number of degrees of freedom (dof) added to the standard 
adiabatic model are shown, f Constraints on a are strongly 
dependent on the chosen pivot point in this case, whereas r^so 
is pivot independent, and therefore provides a good measure 
of the isocurvature contribution. The a limit is shown here 
solely for completeness. 



Repeating the analysis with the Zy parameters we find 
95% upper limits on the isocurvature fraction in terms of 
CMB power, r.so, of 0.13 (CI), 0.08 (NID) and 0.14 (NIV) 
compared to 0.23, 0.13 and 0.24 for the first year WMAP 
data [i^. These constraints from 3 years of WMAP 
therefore show a marked improvement, being ~ 60% of 
those obtained with the first year WMAP data; the im- 
proved polarization data prefer a lower level of isocurva- 
ture. 

The data is fully consistent with a = 0/riso — 0, with 
the goodness of fit — 21n>C improved by only ~ 1 for 
each case. These additional degrees of freedom, however, 
cause the baryon density and spectral index mean val- 
ues to move more than la from their adiabatic values: 
both values are increased by 2a when the NIV mode is 
included, exploiting the degeneracy observed in (43| . 

Our investigation finds that the results are sensitive to 
the choice of prior: constraints on a obtained by sampling 
a and (3 directly differ from those derived from the dis- 
tribution sampled using the parameterization. This 
is demonstrated in the top row of Figure [H where the 
two methods are compared. We see that there is more 
phase space available for models with larger a when sam- 
pling with a uniform prior on the observable isocurvature 
CMB power, than there is when sampling with a uniform 
prior on a. The likelihood of the best-fitting models are 
not affected by the choice of prior however, and we can 
expect the dependence to be reduced as data improves. 
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FIG. 2: 68 % and 95 % 2-dimensional constraints on the am- 
plitudes of generally correlated isocurvature modes, for the 
CI mode (left), the NID mode (center), and the NIV mode 
(right) for WMAP plus SDSS and SNLS data. The top pan- 
els show primordial amplitude contributions in terms of a, (5, 
using flat priors on Zij (line contours) and a, /3 (filled con- 
tours). The {q, 2/3-\/q(I — a)} parameter space is contained 
within a circle of unit radius shown by the dashed line. The 
lower panels show the observable CMB power contributions 
in terms of Zi-j. 
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FIG. 3: The effect of varying the CDM isocurvature spectral 
index independently for an AD -I- CI scenario with WMAP plus 
SDSS and SNLS data: 68 % and 95 % constraints on the 
AD-CI cross-correlation /3 (left panel), and adiabatic spectral 
index (right panel). The dotted line in the right panel shows 



a large phase space is opened up for models with large 
n.iso, within the range allowed by the prior. The 1-D 
marginalized constraints on the isocurvature contribu- 
tion are increased from nso < 0.13 to < 0.15 (95% 
C.I.). The isocurvature tilt cannot be constrained by 
WMAP-l-SDSS-hSNla datasets alone, however ^ show 
that additional Ly-a data prefer higher tilts of 1.9 ± 1. 
When spectral indices are able to vary freely, Viso is a 
good measure of isocurvature because a then becomes ex- 
tremely sensitive to the pivot point at which the spectral 
indices are defined. The scenario we investigate here with 
CDM isocurvature and Uiso riadi is a case in point. For 
models with high niso{0.05/Mpc) the isocurvature power 
on larger cosmological scales is significantly reduced for 
a given a and therefore a is able to be increased to com- 
pensate. The relative power in isocurvature, however. 



rr- 6000 
3. 

4000 

Csl 

2000 



AD 
CI 

AD.CI 



■r, =3% 



10" - 



10^ 



Csl 



cj 10' 
10"^ 



TT — N^.^ 


. TE 






EE ^-—^ 














^ Adiabatic \ 








10 


100 


1000 



FIG. 4: AD+CI: (Top) CMB temperature power spectrum 
for the best-fit model with correlated AD -I- CI modes and in- 
dependent spectral indices (solid) for WMAP plus SDSS and 
SNLS data. The model has nso = 3.4%, a = 0.49, /? = -0.26, 
nadi ~ 0.97, Uiso ~ 2.7, and has —2 In jC lower than the best-fit 
adiabatic model. The contribution of each mode correlation 
to the total spectrum is shown, including the WMAP data. 
Bottom: the CMB spectra are compared to the pure adia- 
batic best-fit model (dashed) . Only the £ < 500 section of the 
TE spectrum is shown. 



roughly indicated by ~ a/ {I — a) x (k/ko)""" is 
not increased. 

The best-fit model is shown in Figure [H has Uadi = 
0.97, Uiso = 2.7 and may be distinguished from the adia- 
batic model at large scales in polarization. The goodness 
of fit is improved, compared to the adiabatic model, by 
a A(— 21n£) of 4 to 11379 for 2 additional degrees of 
freedom. Such an improvement is driven by the use of 
the 3-year WMAP data in the analysis. Future measure- 
ments will help distinguish these, currently degenerate, 
high isocurvature tilt models. 

This subset of models may inform us about double in- 
fiation scenarios. If there are multiple fields driving in- 
fiation, it is possible to generate entropy perturbations 
in which the additional light fields modify the curvature 
perturbations on horizon scales, and also modify the con- 
sistency relations relating scalar to tensor modes [H^, ■ 
In the most general case Uadi niso 7^ ricorr, with each 
related to the slow roll parameters along the flat direc- 
tions of the scalar fields. For example, for theories in 
which one scalar field plays a dominant role in driving 
inflation, but two fields play a significant role during re- 
heating, one finds that riiso ~ ricorr [53j . Analyses in 
which Ucorr has been allowed to vary have found it is 
a nuisance parameter unconstrained by data [3§|, so we 
believe it is reasonable to assume that fixing the scale de- 
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FIG. 5: AD+CI+NID: 2-dimensional constraints show the 
degeneracy between isocurvature cross-correlation amplitudes 
Zij, and the adiabatic amplitude zaa, that allows the destruc- 
tive interference of isocurvature spectra. 



pendence of the cross spectrum , ricorr = {nadi + niso)/'^, 
does not unduly bias the conclusion. For two field infla- 
tion of two minimally coupled scalar fields of mass m/i 
and mi, the magnitude and correlation of the resulting 
CDM isocurvature component are dependent on the ra- 
tio of the masses R = nih /mi and number of e- foldings 
Sk- A bound on R comes from the magnitude of the 
cross-correlated spectrum fs^ 



2/3Va(l-a)|, 



1) 



(i?2 



1 



(11) 



Assuming Sk = 60 we find an upper bound on the ratio 
of the two scalar fields of R < 3.5. This is weaker than 
the i? < 3 at 95% c.l. obtained with the inclusion of 
Lyman-Qf data [4^ , although we caution that constraints 
on a in this case are strongly dependent on the choice of 
prior and pivot scale. 

It would be interesting, but beyond the scope of this 
paper, to place constraints on specific double inflation 
models in which model-dependent predictions for each 
mode's spectral index are included. 



V. GENERALLY CORRELATED 
ISOCURVATURE: MULTIPLE MODES 

In this section we consider models with two additional 
correlated isocurvature modes, (CI-fNID, CI-I-NIV, 
NID-I-NIV), and finally a model with the full set of adi- 
abatic and three isocurvature modes. We sample the 
modes using the parameterization given in equation ([7]) 
for N — 3 and N = 4. Table [11] shows constraints 
for the relative mode contributions for this set of mod- 
els. We also give the primordial amplitudes An of the 
auto-correlations contributing to the best-fitting models, 
where Ce = J2fj=i AijC}^ . 

Two isocurvature modes: As is shown in Table HIl mod- 
els with two modes permit far more isocurvature than 
those with a single mode. 

Although when considered individually, the neutrino 
velocity isocurvature modes allow the largest isocurva- 
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FIG. 6: AD+CI+NID: Spectra as in Figure H The model 
shown has correlated modes, with rtso ~ 41%, and cosmo- 
logical parameters Q,i,h^ — 0.026, Qch^ = 0.12, Qa = 0.73, 
T — 0.10, Ha = 0.92, bsDSS ~ 0.98. The isocurvature spectra 
add destructively, canceling almost completely. 



ture fraction, interestingly when two modes are included, 
joint CDM and neutrino density isocurvature (CI-I-NID) 
allow the most freedom, more than twice as much isocur- 
vature {riso = 0.4±0.1) as the combinations including the 
neutrino velocity mode. This freedom arises from degen- 
eracies within the isocurvature cornponents destructively 
interfering, originally observed in [40| and shown in Fig- 
ure [5] Degeneracies with the NIV mode do also exist 
where the spectra add constructively, but such models 
have large baryon densities ruled out by current BBN 
measurements. 

Figure [5] shows the best-fit CMB temperature spec- 
trum for the CDM + neutrino density isocurvature model 
with risa = 41%, compared with 47% with WMAP year 
1 data 40]. The contributions from all six correlations 
(three auto-correlations and three cross-correlations) lead 
to greater large scale polarization CMB power, but cancel 
almost completely in the CMB temperature and galaxy 
power spectrum. 

All three two-mode models prefer baryon densities 
higher than the concordance value (with mean values 
0.025 < rjfe/i^ <- 0.027), despite the BBN constraint. 
The spectral index is also more poorly constrained, 
with the CI-I-NID models preferring a low spectral in- 
dex (0.93 ± 0.03), and the CI-ENIV (0.99 ± 0.04) and 
NID-I-NIV (1.03 ± 0.03) preferring larger values. The 
other cosmological parameters are consistent with adia- 
batic A CDM. 

Three isocurvature modes: When all three modes are 
included, the constraints are highly sensitive to the BBN 
constraint due to the strong degeneracy between fl^h'^ 
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FIG. 7: AD+CI+NID+NIV: Spectra as in Figure H The 
model shown has Viso ~ 62%. No BBN or bias constraint has 
been imposed; the cosmological parameters are Qbh'^ = 0.037, 
Qch'^ = 0.13, riA = 0.75, r = 0.10, n, = 0.98, bsDSS = 1.2. 



and the NIV amplitude. With the BBN prior and SDSS 
bias we find nso = 0.44±0.09, increasing to nso = 0.51± 
0.09 when no BBN or SDSS bias priors are included, to 
be compared to 0.57±0.09 found for the first year WMAP 
data |43l |. 

A well-fitting model (A(— 21n£)=-6 with 10 extra de- 
grees of freedom in comparison to the adiabatic best fit) 
with a majority of the power coming from isocurvature 
modes, riso = 62%, is shown in Fig. \7\ This model was 
obtained without the BBN prior and has a large baryon 
density {flth'^ — 0.037); however even including the BBN 
constraint, the baryon density is higher than the con- 
cordance value for this class of models (0.031 ± 0.003), 
and Qch'^ is raised by Icr to 0.124 ± 0.007. Figures M 
and[7]indicate that precision small-scale temperature and 
large-scale polarization will to more tightly determine the 
underlying initial conditions. In particular future small- 
scale CMB experiments should strengthen the constraint 
on the baryon density, since high baryonic isocurvature 
models are more degenerate at smaller scales than their 
low baryon counterparts. 

Given that the data can still only poorly constrain 
the isocurvature contribution for these multi-mode mod- 
els, the constraints presented here depend on the prior 
distribution we have chosen. Previous work has shown 
that other parameterizations can decrease or increase this 
contribution [4^, depending on the phase-space volume 
available for purely adiabatic models compared to mixed 
isocurvature models. As with the single mode case, im- 
proved data will help limit this prior dependency. 



VI. CONCLUSION 

We have investigated the constraints on the presence of 
a variety of isocurvature modes in the initial conditions 
of structure formation, in light of recent observations of 
temperature and polarization CMB data and large scale 
structure and supernovae surveys. 

The improved WMAP data, with the inclusion of low I 
polarization measurements, has strengthened these con- 
straints on the contributions of individual isocurvature 
modes, with the polarization data disfavoring models 
with a large isocurvature fraction. Scenarios with ei- 
ther CDM (or baryon) or neutrino isocurvature allow 
only a very limited contribution, which can be translated 
into strong constraints on the curvaton model and some 
double-field infiationary models. Although models with 
multiple isocurvature modes do not offer a significantly 
better fit to the data, models with non-zero isocurvature 
fluctuations fit the data as well as the adiabatic model 
and can comprise the majority of power when additional 
modes are considered simultaneously. 

Of the models with more than one isocurvature mode, 
those most likely to pose the greatest difficulty for distin- 
guishing with future data are those with large fractions 
of both correlated CDM and neutrino density isocurva- 
ture, which provide the best fit to the data, and due to 
their destructive interference are highly degenerate in the 
CMB and galaxy power spectra. Those with neutrino ve- 
locity fluctuations (both two and three mode models) are 
better constrained by BBN and bias measurements. 

With WMAP plus LSS and SN data, the baryon den- 
sity and spectral tilt are found to be sensitive to the 
inclusion of isocurvature modes. With the current data, 
the reionization optical depth however is robust despite 
the modifications isocurvature models can make to large 
scale polarization spectra. Extending beyond the A CDM 
scenario, given results found in [4^, we would not expect 
that allowing independent tilts for all the modes, or in- 
cluding curved geometries, to have a large effect. Future 
B-mode polarization data will help break degeneracies 
between tensor and isocurvature modes that would cur- 
rently arise from large scale temperature CMB data. 

Future small-scale temperature and polarization data, 
together with improved galaxy and Lyman-a power spec- 
trum measurements, should help constrain a subset of 
the models we have considered, but improved large scale 
CMB polarization data from WMAP and in particular 
Planck, demonstrated in [s^l , will be crucial if we are to 
strongly constrain this general set of correlated isocurva- 
ture models. 
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